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Abstract 

We report here our lattice simulation on the charmonium spectra in the 
quenched approximation. Because the full adjustment on the nonperturba- 
tive parameters such as Ce, Cb, mo a s and r s needs many calculation time, we 
only adjust two of them, m a s and r s ( £3 plays this role in the paper) but with 
some rescale for mass splitting. After the rescale, we find that our results are 
in agreement with the experiment ones. 



1 Introduction 



Due to the importance of the charmonium spectra in particle physics, there are 
many jobs on it in lattice QCD. These jobs did their ways with NRQCD approach^ 
or with Wilson action with clover-improvement Pi PI [4] [5] _ 

We report here our lattice results on the charmonium spectra on anisotropic 
lattice ( we will also present our results on the bottomonium spectra in the future). 
Paralleling with references [2] -[5], our results were obtained with the Sheikholeslami- 
Wohlert action® '^(SW action) as well as with the tadpole improvement. As many 
authors pointed out, such improvements make us get interesting results with the 
fewer calculation. In fact, our simulation was done on personal computations. 

As pointed out in reference [13], if we adjust parameters to take the right hy- 
perfine splitting in this action, the mass splitting generated by the orbit excita- 
tion, including orbit angular momentum and therefore spin orbit interaction, will 
be incorrect (we call this mass splitting s.o. splitting thereinafter). Therefore, to 
get the correct charmonium spectra, or the mass splitting, some argument on the 
mass dependent choice of C E was done in the article [13]. In fact, the choice of the 
clover coefficient have been performed in various ways, which include applying differ- 
ent orders in perturbation theory^ and introducing non-perturbative improvement 
criteria using Schrodinger function®. Those schemes, however, need much more 
computer or calculation time. 

If parameters on lattice are not maladjusted so much, then the mass splitting 
and the gross mass, which is mainly dictated by the static quark mass M 1 , are in 
different scale. To make the scale at the same, maybe some rescale scheme for the 
measured mass splitting is also proper. After the rescale, we find that our results 
are in good agreement with the experiment ones. 

The section 2 shows some notations we should considered when we begin to 
calculate, and the section 3 and 4 is our simulation detail and results. At last, we 
present the summary in the section 5. 



2 The Choice of Parameters for Simulation 

To perform the simulation, we adopt the tadpole improvement action for gluon as 
follows: 



c _ a(5Psp . 4^P tp 1 P sr 1 £P s tri /, x 



where P sp {P s t) is the sum of all spatial(temporal) plaquettes, P sr represents 2x1 
spatial rectangular Wilson loops and P str is the Wilson loops with 2 spatial and one 
temporal links. £ is the aspect ratio: 



This gluonic action, as many authors pointed out^, may improve the action up to 
0(aj,a 2 t ). 

For quarks, we adopt the SW action. This action was proposed initially to 
improve the continuum behavior for the fermions. But as suggested by liischer and 
Weiszt 11 !, one had better take the scheme with the on-shell improvement. This 
scheme had been illuminated systematically in reference [13] in action form as well 
as in Hamilton form. Here we shall use the scheme to determine some parameters 
in SW action. 

We take the SW action as: 

S g = E $( x ) W* s + \ [674 V 4 + E 7iVi(l -\ A,)- ^1 A 4 

~ E( A * - A a?) + CeJ2 + c B y E 

V^{x) = u^{x)^{x + fi)-ul(x-fi)4)(x-fi), 

A^(x) = u^x^ix + ft) + ul(x - p)ip(x - p) -2ij}{x), (3) 

where mo is the bare mass of the quark. We shall determine the parameters £i, £3, 
Ce and Cb in Eq. (3) in the following. 



2.1 The choice of parameters £1 and £3 

The propagator for quark is 

<^(t',p')^(t,p) >= (27t) 3 5(p'-p)C(t>-t,p), (4) 

where 

with the applying of the residue theorem. In Eq. 5, Z(p) is the function which is 
independent of time t. The energy E of the quark with momentum p is determined 
by 

G 2 + w 2 + ^ = 2G£i cosh Ea t , 

G = m + ^1 + — E(! - cos ^ + -(1 - cos p^ 2 ), 



w 



1 E( sin ^ ( 4 - cos ^)) 2 - ( 6 ) 



9 ^7 



To get above equation, we set a s = 1. So, one may expand G and u> 2 according to 
the power of pf 

2 

G = m + 6+E^ + °(P 6 )' 

= Y.p*+o(p 6 )- (7) 



At mo = 0, one gets: 



(bEa t ) 2 + 0(at) = J2Pi+°(ptl (8) 

i 

where at = l/£. Because of the Eq. (8), we choose £i = £ with the tuning of the 
parameter £3 to make the following dispersion relation of the particle ifj be satisfied: 

E 2 (J/4,) = m 2 (J/4,) + c 2 Y,Pl (9) 

i 

where c = 1 is the effective velocity of photon(EVP). 

One can not make all the dispersion relations of the calculated particle be sat- 
isfied at the same time with the tuning of the sole parameter. Phenomenally, this 
infraction manifests the fact that the relative movement of two constituent quarks 
in meson is various with different meson( See appendix for more detail). 



2.2 The Choice of C E and C B 

The key points of the choice for the clover coefficient C E is the isospectrum transfor- 
mation and the redundant coupling. After some algebra calculation, one can easily 
see that the choice C E = 1 is proper^. 

There are a few choices of C E . For instance, one may choose C E = 1 or C E = 
[12] i n NRQCD. In fact, following the method developed in reference [13] and [3], 
one may see that in this action the parameters should satisfy 



C B = 1, 



£ / 1 2 w 



A = — 5 — < 1 _ ) 

l+m t y \n(l + m t ) m t (2 + m t ) J 

C - j I ( £ 2 m t {2 + m t ) 

where m t = The above condition is at tree level but with tadpole improvement. 

While for getting the nonperturbative parameters we should adjust these parameters 
with the same physical principle that led to Eq. (10). For instance, to get the 
nonperturbative £3 we adjust £3 to make Eq. (9) to be satisfied. 

We select Cb = 1 for our simulation. While noticing that C E = £/^3 at vanished 
m t we just take a mass independent setting for C E for simplify: 

C F = i. (11) 



6 



3 



The choice in Eq. (10) and (11), however, is at tree level and sometimes a little 
maladjusted. To compensate the maladjustment will lead to a mass splitting rescale. 

Thus we have only two parameters £3 and mo to be determined in the simulation: 
one is adjusted to give the correct energy-momentum relation for J/ip and the other 
will be rectified to determine the lattice spacing, a s (a t ). 



3 The Detail for the simulation 



In this work, we calculate five meson masses of charmonium, named as l So(r) c ), 
3 Si(J/ip), 1 Pi(h c ), 3 Po(Xco) an d 3 -Pi(Xci) and their excited states. The choice of the 
operator for definite meson is as the standard one' These operators are called T 
ones in ref. [3]. 

To determine £ 3 , we calculate the energy for four low-lying momentum j^pa s = 
(0, 0, 0), (0, 0, 1), (0, 1, 1) and (1, 1, 1) of the particle J ftp and fit them by 



r 2 



(E(\p\)a t ) = (E(0)a t y+^pp\a s y. (12) 

We tune £3 to make the relation Co = 1 be satisfied in one percent accuracy. 

There are types of scheme to determine m ,or equivalent a t . For instance, one 
may determine m based on the hyperfine splitting or based on the s.o. splitting. 
As for the first type, One may take the ratio of the masses of S-state mesons: 

m{ib')m(r) c ) _ . . 

0.987. (13) 



m(r)' c )m(J / ip) 



And as for the second, one may take, for example, the ratio of the mass of Xci and 
J 1 4 

m(Xcl) - 1.134. (14) 



m(J/ijj) 



Then, we can extract a t by 



= m(lS)at 
1 m(lS) eX p 

where m(lS) = (m^S) + 3m(l 3 S))/4 is the measured average mass of l 1 ^ and 
1 3 S state on the lattice and m(lS) exp is the experimental one, 3.0676Get> . 

In this paper we iterate mo and £3 to make both the equation (12) and (13) be 
hold at the same time. 

Because of our computational limitation, we did our simulation on a L x T = 
6 3 x 36 lattice with £ = 6.0. Parameters for the simulation are displayed in table 1. 

For a definite T operator we measure its correlation function as: 

c state (t, P ) = E < ^(x, t)rv(x, t)$(o, o)r^(o, o) > e ^ p - x . (ie) 

X 

The energy can be extracted by n-hyperboliccosine ansatz( n is the number of the 
extracted masses and T is the periods in the time direction): 

n— 1 

C state (t) = £ a,(e""^ + e~ m ^), (17) 
3=0 



which means that we should find parameters a(j) and m(j) to minimize the function 



If n—l 

£ Wi(C(ti) - £ aj (e- m * + e- m ^ T -^)) 2 /C(t t ) 2 , (18) 

i=h 3=0 

where W{ oc a is the weight of the correlation C(ti). 

We choose n = 2 for all the measured particles in the work. Apparently, the 
excited states we measured here is not the real first excited states in nature but the 
mixture, the mainly component of which is the first and the second excited states. 
Thus all the measured excited masses are larger than the experiment ones. We 
expect that the scheme we adopted in Eq. (14) would decrease this mixing impact. 

For each U and If in Eq. (18) one may find a group of corresponding parameters 
rrij and dj. Then, one may choose a set of parameters k, If, rrij and a,j, which make 

the function XX - T" "I t) § e ^ ^ s minimal points for all the sets. But due to the 

fact that nearly all the minimal point for a definite operator is approached at the 
set in which k — 1 and If — T — 1, we fix our set in which k = 1 and If — T — 1. 
This is as expected, for the more measured data means the more information if we 
excluded the case of U = or If = T, in which the counterterms will enter in. 



4 Simulation Results 

In table 2, we list the results of the typical dispersion relation and the fitting data 
for low-lying particles at (3 = 3.0 and m a s = 1.06. 

From the simulation we know that the EVP for the particle r] c is equal to the one 
of the particle J/ip in about 2 percent. This is obvious from the appendix since their 
mass splitting is mainly due to the spin splitting^. The similar result is obtained 
by other authors^'~^. 

On the other hand, we find that the EVP in different mesons usually have the 
relation: 

c(h c ) < c(xci) < c(xco) < c(J/V>) ~ c(rj c ). (19) 

This possibly implies that the expanding point p s, or roughly saying, the main 
distribution of the relative momentum, can be estimated as 

Po(K) > Po(Xci) > PoiXco) > PoWO ~ Po(Vc)- (20) 

From the data we see that the EVPs for all particles except meson r/ c and J/t/j 
are usually smaller than unity, which means that the kinetic mass M 2 is larger than 
the static mass M\{ see detail in ref. [13]). At the same time, to get the correct 
Xj mass splitting, ref. [13] suggested a mass dependence of Ce, which leads to a 
scheme of adjustment of three parameters: m a s , £ 3 and Ce- 



To combine the scale of Mi, M 2 and Me and decrease the number of adjusted 
parameters, we adopt such progress: we adjust only two parameters m a s and £ 3 
but with the mass splitting rescale. 

We take Eq. (13) to determine m a s and measure the mass of particle i] c and 
J ftp, but for other particle, since we should consider the s.o. splitting, we redefine 
its mass as follows. 

We rescale the mass splitting between 1 3 P X and averaged IS* state as 

m(xd) ~ m(lS) = O.U3Gev. (21) 
Therefore, the modified mass for a definite meson is 

m mod = 0.443 p -^(^H + m{1§) {Geo) (22) 
m{Xci)a t - m(lS)a t 

where ma t is the dimensionless measured mass on the lattice and a t is extracted by 
Eq. (15). 

The simulation parameters are listed in table 3. Our results for scale are not 
inconsistent with the result obtained from the Sommer scale r in the calculation of 
the gluon spectra^ although they are achieved in very different aspect. 

We list our modified charmonium spectra in table 4, which is, as well as in other 
references, in agreement with the experiment. From the data we find that at the 
finite lattice spacing the ratio is smaller than the experiment one. To get 

the correct ratio one may decrease m and therefore achieved a decreased lattice 
spacing, a s . An estimation for the decreased itlq shows that it approximately leads 
to the same result as in reference [17]. To decrease the computer time, we take the 
scale determined by Eq. (13) with the notice that the ratio between a s obtained by 
us and that obtained by C.J. Morningstar is nearly the same at different f3, i.e., the 
ratio is 1.6 at (5 — 3.0 and 1.7 at f3 — 2.5. 

While compared with our results with those in [5], one may find, that our result 
for hyperfine splitting is smaller than that in [5]. Maybe the smallness is caused by 
the too large lattice spacing and/or the contamination of hyperfine splitting by the 
s.o. terms in the determination of m a s in Eq. (13). 

For the known excited states, our results, as those obtained in ref. [3], are larger 
than the experiment ones. But for unknown states, it needs some more explore to 
explain why our result are also larger than those in ref. [3]. 

While following the reference [13], one should do this rescale between xj states. 
But, from the results we find that our rescale scheme also obtains the correct xj 
states. 



5 Conclusion 



For the simulation of the charmonium spectra, because the ratio between the rest 
mass Mi and the kinetic M 2 are not always the same, so the scale for mass splitting 



Am and the the gross mass m is different. Since Am is not mainly dictated by the 
static mass M 1: reference [13] suggested the following strategy: forget about Mi and 
adjust the bare mass so that the kinetic mass M 2 takes the physical value. 

To make mass splitting and the gross mass at the same scale, one may adjust 
parameters or rescale mass splitting. This rescale combines the scale of the gross 
M and the mass splitting, AM. After the rescaling, we get the results which are in 
agreement with the experiment ones despite a little maladjustment of the parame- 
ters. 



Appendix A 



From phenomena, one may regard the meson consist of one couple of constituent 
quark and antiquark and one can also in principle write the wave function for (anti- 
) charm in momentum and spin space. The relative movement between them depicts 
the different mesons. For instance, the orbital angular momentum in r\ c is S-wave 
while that in h c is P-wave. On the other hand, for cc system the typical velocity in 
particle is v 2 ~ 0.3, i.e. the energy scales m c v 2 /2 is about 200-800 Mev. This implies 
that the mass splitting generated by the relative movement between quark and 
antiquark is mainly dictated by the dynamical behaviour of quarks at the vicinity 
around m c v 2 /2. So, maybe the better expanding points in meson is p 2 ] (/2m c ) rather 
than pi = 0, where p\ is phenomenal parameter to describe the relative movement 
of quarks in mesnon. 

In the following we only consider the effect of the relative movement in 1+1 
dimension lattice. 

In continuum, the energy-momentum relation for the quark is: 



E 2 (p) = m 2 + ^f- = m 2 + ^(p 2 + Ap 2 ), (i) 



where all the quantity is dimensionless one, i.e. E is Ea t , quark mass m is ma t and 
p is pa s . One may rewrite it as: 



e E + e -E = e ma + e ~ma + ^P^ma _ e ~ma^ + QUAp 2 ) 2 ), (fi) 

2ma^ 2 



„2„2 



where a 2 = 1 + = 1 + l 2 v 2 . One may redefine the phenomenal parameter v 2 
to set 7 = 1. 

On the other hand, on an 1+1 dimension lattice we have: 
w 2 = wl + ^Ap 2 , 

2 ^ • 2 I a \2 

w = -an p (4-cospo) , 

sinp (4-cospo), /1 v 

fi = ^ (4cosp - cos2p ); (m) 

9p 



and 



G = Go + f 2 Ap 2 + 0((p 2 y), 

2 1 
G = m + £+^sin 2 po/2(l + -sin 2 p /2), 

sinpo,, v ( . x 

h = — (5-2cosp )- (iv) 



Therefore, according to Eq. (6), on lattice one gets 

e E + e~ E = g° + 1- + # + £ + A - + 4))^P 2 + 0((Ap 2 ) 2 ). (v) 



Comparing with Eq. (ii)and setting c = 1, one knows that 

c = l-r 2 (l + ti;g/f) 

S3 sinh met /i J j 



where 



ma t 

2 

= Z(t-1-— g), (vi) 

«3 



# = sin 2 p /2(l + ^sin 2 p /2), 



&+v / & 2_ 4(1+w 2 /a 
* ~ 2 

b = e ma + e- ma . (vii) 

While noticing that functions / 2 , /i, Go an d w o are even ones f° r Po, one ma y 
coarsely replace p with in those functions. We expand Eq. (vi) with respect 
to ma s : 

m a s = ma s + 0((ma s ) 3 ) (viii) 

To applying the above equations, we should consider the intricate additive and 
multiplicative renormalizationt-'-'^ . 

The above argument is somewhat pedagogical, but it manifests the correct qual- 
itative behavior, even on the 3+1 dimension lattice (see detail in the paper). For 
instance, from Eq. (viii) it is easily understood that the error of the pole mass will 
be proportion to ma if we choose £ 3 = £ naively. 
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p 


< 


sweep/confg. 


confgs. 


2.5 


0.423 


40 


100 


2.8 


0.463 


40 


100 


3.0 


0.486 


40 


100 



Table. 1 Parameters for the simulation with the scale set by Eq. (13). 



(£IpK) 2 





1 


2. 


3 


Vc 


0.276(3) 


0.309(2) 


0.340(2) 


0.369(2) 




0.278 


0.308 


0.339 


0.370 




0.286(1) 


0.317(1) 


0.349(1) 


0.3775(13) 




0.286 


0.317 


0.348 


0.3786 


h c 


0.409(13) 


0.400(10) 


0.416(8) 


0.432(7) 




0.398 


0.408 


0.419 


0.429 


Xcl 


0.401(9) 


0.414(8) 


0.432(7) 


0.448(7) 




0.400 


0.416 


0.432 


0.448 


XcO 


0.380(11) 


0.405(11) 


0.428(12) 


0.450(12) 




0.380 


0.404 


0.427 


0.451 



Table. 2 {Ea t f v.s. (^) 2 for different mesons. The data in upper row are measured 
ones and the data in lower row are fitted ones. The fitting parameters a and b of the 
formula (Eat) 2 = a + b ^ p ^ for different meson are 0.2776(20) and 1.007(34) for r] c , 
0.2863(12) and 1.010(17) for J/ip, 0.398(7) and 0.34(11) for h c , 0.3995(58) and 0.530(93) 
for Xcl, and 0.380(8) and 0.77(13) for XcO respectively. The corresponding velocities of 
photon for different meson are: 1.003(17) for r/ c , 1.005(8) for J/if>, 0.58(10) for h c , 
0.73(6) for Xcl and 0.88(8) for XcO respectively. 



p 


& 


m a s 


a s (fm) 


La s (fm) 


EVP(J/^) 


EVP(r ?c ) 


2.5 


3.31 


2.33 


0.322 


1.93 


1.000(10) 


1.022(21) 


2.8 


3.92 


1.57 


0.258 


1.55 


1.006(7) 


1.018(15) 


3.0 


6.38 


1.06 


0.208 


1.25 


1.005(8) 


1.003(17) 



Table. 3 Parameters for the simulation with the scale set by Eq. (13). 



state 


f3 = 2.5 


(3 = 2.8 


f3 = 3.0 


a s -> 


Exp.(Gev) 


Vc 


3.030(11) 


3.032(11) 


3.030(14) 


3.030(22) 


2.980 




3.080(47) 


3.0795(44) 


3.080(7) 


3.080(10) 


3.097 


h c 


3.553(135) 


3.500(68) 


3.539(43) 


3.546(117) 


3.526 


Xcl 


3.511(95) 


3.511(50) 


3.511(32) 




3.511 


XcO 


3.471(139) 


3.447(73) 


3.436(46) 


3.412(119) 


3.415 




4.309(26) 


4.347(17) 


3.982(15) 


3.766(31) 


3.594 




4.316(11) 


4.359(7) 


3.993(6) 


3.795(13) 


3.686 


K 


4.748(107) 


4.964(68) 


4.633(58) 


4.564(126) 




Xcl 


4.889(92) 


5.079(57) 


4.689(45) 


4.534(102) 




XcO 


4.980(150) 


5.145(94) 


4.721(70) 


4.510(161) 





Table 4 Results of the charmonium mass M and the mass splitting AM 



in units of Gev at £ = 6.0 with the same scale setting as in Table 1. 



